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Abstract. In order to study the influence of the environment 
on the stellar population of early-type galaxies we have an- 
alyzed the relations between the density of the environment 
log(p) and the residuals to the Mg 2 - Co relation, R a , and be- 
tween log(p) and the residuals to the Fundamental Plane, Rf. 
Our sample of galaxies covers the range of densities, between 
field galaxies and poor clusters. We confirm the existence of 
a slight environmental bias on the mean stellar population and 
we show that this effect is due to a small fraction of galaxies 
containing a younger stellar sub-population while the majority 
of galaxies, in any environment, are equally old. The galax- 
ies which are likely to contain a young stellar sub-population 
(negative R a and Rf) are preferentially found in the less dense 
environment. 

This segregation of the stellar population may reflect the 
morphology-density relation or may result from a higher rate of 
recent star formation activity in sparse environments. Using the 
rotational support, V max /ao , to assess the presence of a disk 
and to parameterize the morphology, we subtracted the contri- 
bution of the morphological segregation to population segre- 
gation. It is not yet possible to rule out that the morphology- 
density relation is at the origin of the population-density rela- 
tion. 

Key words: galaxies: general - galaxies: elliptical and lentic- 
ular, cD - galaxies: fundamental parameters - galaxies: stellar 
content - galaxies: formation 



1. Introduction 

The comparison o f the s tellar populations at diff erent distances 
(e.g., Bender et al. |1993|, Ziegler & Bender |1997|, Stanford et al. 



1998) shows that early-type galaxies (E and SO) are essentially 



old and undergo only a passive evolution. Even in sparse envi- 
ronments, where hierarchical models of galaxy formation pre- 
dict a younger age, the mean ages are not dramatically younger 



Send offprint requests to: Ph. Prugniel (prugniel@galaxies.univ- 
lyonl.fr 
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(J0rgensen 1997 Bernardi et al. 1998). However, a fraction of 



elliptical galaxies contains also a younger component. This is 



apparent from statistical analyses, eg. Forbes et al. (1998), and 



strong evidence comes from the observation of merger rem- 



nants (e.g., NGC 7252, Schweizer 1982) which will turn into 



"normal" ellipticals after a couple of gigayears. 

These intermediate-age sub-populations provide the frame- 
work for understanding the connection between^zne structures 



and spectral peculiarities (Schweizer et al. 1990) and probably 



also for interpreting the peculiar velocities found in clusters 



(Gregg 1992) and residuals to the Fundamental Plane (Prug- 
niel & Simien[l99§). 

It is generally accepted that the rate of star formation in 
early-type galaxies is enhanced in low-density environment 



(Schweizer & Seitzer 1992, de Carvalho & Djorgovski 1992, 
Guzman et al. [1992 Rose et al. 1994, J0rgensen & J0nch- 



S0rensen 1998| ). It results in a population - density relation: 
The metallic features in the spectra (eg. Mg 2 ) are weaker and 
the Balmer lines stronger in low-density regions. 

However, the interpretation of this environmental depen- 
dence is not straightforward. For example, Mehlert et al. ( 1998[ ) 
interpret the dependence of Mg 2 on the distance to the center 
of the Coma cluster, previously found by Guzman et al. (1992), 
as a bias in the morphological classification. The "young" 
early-type galaxies found in the outskirts of Coma would ac- 
tually be lenticular galaxies, and the population segregation 
would simply reflect the morphological segregation. 

Because the stellar population varies with the Hubble type, 
at least part of the population - density relation is likely a 



by-product of the morphology-density relation (Dressier 1980 



1984, Whitmore et al. 1993; the fraction of early-type galax- 



ies is higher in dense regions). Conversely, the population seg- 
regation may bias the morphological classification (for spiral 
galaxies) and hence contri bute to the morphological segrega- 
tion (Koopmann & Kenney 1998 ). 

The question left open is to determine the fraction of the 
population segregation due the morphological segregation. It is 
connected with that of the origin of the morphological segrega- 



tion (Martel et al. 1998). This segregation may be due to initial 



conditions, morphological evolution (due to gas stripping and 
mergers), or a combination of both. The first hypothesis is dif- 
ficult to defend in the light of the existence of clear cases of 
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merger remnants. If the environment is responsible for a sig- 
nificant morphological evolution, it will also be at the origin 
of a population segregation (at a given morphological type) be- 
cause the morphological transformations are accompanied with 
star formation. 

In elliptical galaxies, the "young" stellar sub-populations 
are likely due to merging with a gas-rich companion that oc- 
curred at maximum a couple of gigayears ago. In lenticular 
galaxies, the recent populations may be the result of the resid- 
ual star formation in the gaseous disk. The first class of ob- 
ject will not be subject to the morphological segregation, while 
the second will be (in this case the morphology is related to 
the importance of the disk). Subtracting the contribution of the 
morphological segregation to the observed population segrega- 
tion, would in principle allow to determine the present rate of 
environmentally triggered star formation. 

The aim of this Letter is to study the population segrega- 
tion in the sample of nearby early-type galaxies presented in 



Prugniel & Simien (1996), and to address the question of its 



origin. We will first parameterize the density of the environ- 
ment using the HYPERCAT database (http://www-obs.univ- 
lyonl.fr/hypercat). Then we will use two diagnostics to detect 
the presence of a young stellar sub-population. The first one 
is the analysis of the residuals, R a , to the well-known relation 
between the magnesium line strength inde x Mg ^, and the cen- 
tral velocity dispersion <7rj (Terlevich et al. 



1981, Bender et al. 



1993). The second approach similarly analyses the residuals to 



the Funda mental Plane (FP) , Rf, f ollowing the line of Prugniel 
& Simien ( \l99j , |l996j and |1997| , collectively referred as PS). 
The analysis of the rotational support of the galaxies where a 
young sub-population is detected will allow to determine if the 
age segregation is linked to the morphology-density relation. 

2. Analysis 

Our sample of galaxies, described in PS, consists in nearby 
early-type galaxies in different environments. Photometric and 
kinematic observations were obtained using the 1.20m and 
1.93m telescopes of Observatoire de Haute Provence and the 
CARELEC long-slit spectrograph. This material is already pre- 
sented in details in Prugniel & Heraudeau (1998), Prugniel & 



Simien (1994), Simien & Prugniel (1997a, 1997b, and 1997c) 



and Golev et al. ( 1999 1. All the data are also available through 



HYPERCAT. Both Mg 2 (Lick system) and gq are aperture cor- 
rected and standardized to an homogeneous system (see Golev 



& Prugniel 1998) 



The sample includes bona-fide elliptical and lenticular 
galaxies as well as merger remnants. The classification was 
made on the basis of literature assessments and presence of 
peculiarities: galaxies with morphological disturbances, dust- 
lanes, post-mergers ... were rejected from the bona-fide ellip- 



tical subsample (Prugniel & Simien 1996). This classification 



is clearly subjective, and some of the bona-fide ellipticals may 
actually hide a disk and should be re-classified as lenticular. 
This is a general drawback of the mo rpholo gical clas sifica tion 
of galaxies (see Kormendy & Bender 1994 Andreon 1998 ). In 



this paper, the bona-fide ellipticals subsample will be used as a 
reference for determining the FP and Mg 2 - <jq relations. 

2.1. Parmeterizing the density 

To study environmental effects, we will define a parameter 
to measure the density of the environment. Ideally, the den- 
sity of the environment should be the number of galaxies per 
megaparsec cubes locally measured in a given volume around 
a galaxy. For our purpose, the smoothing-volume will be the 
group or cluster the galaxy belongs to. The underlying idea is 
to characterize the mean environment of a galaxy over its 10 10 
years life in order to connect the density of this environment 
with the stellar population built over the same period. 

We will start from the sample of galaxies with measured 
radial velocity smaller than 9000 km sec _1 extracted from HY- 
PERCAT (the velocity data compilation primarily comes from 
the LEDA database]]). It contains 22689 galaxies, but is not 
complete in apparent magnitude or diameter. Extracting a 
magnitude-limited whole-sky sample would produce a list of 
galaxies too restricted to allow proper determinations of the 
density of the environment. We prefer to apply completeness 
corrections to the present sample. Since the number of galaxies 
per 1000 km sec _1 shell is almost constant in the sample, as- 
suming a uniform distribution of the galaxies in all the volume 
suggests a completeness correction: p c — p u * (Vdist/1000) 2 
(where p u and p c are respectively the un-corrected and cor- 
rected densities, in relative units, and Vdist the distance of the 
group or cluster in km sec -1 ). We do not apply correction for 
galaxies nearer than Vdist = 1000 kmsec -1 . This correction 
supposes that the list of galaxies with known redshift uniformly 
samples the real distribution of galaxies. This is a crude hy- 
pothesis, since most redshift surveys concentrated on limited 
regions (in particular clusters), and this could lead to an over- 
estimation of the density in the regions of deep redshift sur- 
veys. To check the magnitude of this bias, we compared the 
redshift sample with the UGC and ESO samples (which are 
diameter limited). We found, that the density projected on the 
sky of the redshift catalogue, normalized to the density in the 
UGC and ESO samples, fluctuates by a factor 2 rms (the den- 
sity was smoothed in 5 degree diameter disks). We did not at- 
tempt to use this result to modify the completeness correction, 
because that would have supposed to apply a further correc- 
tion to the UGC/ESO samples, making assumptions on the lu- 
minosity function of galaxies, which would carry uncertainties 
not smaller than a factor 2 also. 

The algo rithm to group the galaxies is described in Golev 
& Prugniel ( 1998 ). It associates a group to all the galaxies in 
the sample, returning the mean (flow) velocity, the radius and 
the number of galaxies grouped. This is used to compute the 
density and also the aperture corrections which are applied to 
the Mg 2 and <tq data. Rescaled to the mean density of Virgo, 
this algorithm gives p = p c /p c (Virgo) = 0.6 for Fornax, 4 for 
Coma and about 0.3 for Leo or NGC 5846 groups. These values 
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restore the hierarchy of concentrations between these different 
groups and clusters, thus validating our measurements of p. The 
density provides a smooth parameterization of the field-group- 
cluster classification. Our sample covers mostly the range of 
low densities, from field galaxies to poor clusters. 

2.2. Residuals to the Mg 2 - Co relation 

The existence of a tight correlation between Mg 2 and ctq was 



discovered by Terlevich et al. (1981) and further discussed in 



Burstein et al. (1988) and Bender et al. (1993). The slope of 



this relation is clearly due to the metallicity, and not age, as as- 
sessed by the constancy of the slop e of th e color-magnit ude re - 
lations out to z«0.9 (Stanford et al. |l998j Kodama et al. |l998| ). 
At variance, the spread around this relation is likely due to the 
contribution of young sub-populatio ns whi ch results in a skew- 
ness of the residuals (Burstein et al. 1988, Prugniel & Simien 



1996) 



We have used the subsample of bona-fide ellipticals to fit 
the relation between Mg 2 and Co taking properly into account 
the errors on both coordinates. The residuals, R a , fitted on 308 
objects, are defined as: 

R a = Mg 2 - 0.225 ± 0.052 x log(cr ) + 0.235 ± 0.083 (1) 

The estimated error on calculated Mg 2 is 0.026. This fit is in 
agreement with Guzman & Lucey (1993) and with Davies et 
al. ( |1987[ ) who both find a slope of 0.20 ± 0.05. 

Separating this sample in two parts according to the den- 
sity, we find < R a >= 0.003 (154 galaxies) and < R a >= 
—0.004 (154 galaxies) respectively for the high- and low- den- 
sity subsamples. The difference, 0.007±0.005, perfectly agrees 
with s imilar estimates by J0rgensen ( 1997 ) and Bernardi et al. 
( |1998| ). 

In the present analysis, we are interested in the nature 
and environment of early-type galaxies hosting a young sub- 
population. For this purpose, Fig. [j] presents the average 
log(p) as a function of R a . The galaxies are grouped in bin 
of similar R a (ie. presumably similar stellar population). The 
trend for galaxies with most negative residuals to be located in 
lower density environment is clear. 

We have over-plotted in Fig. [j] the subset of elliptical galax- 
ies, it is interesting to note that the same tendency is observed, 
but with a slight offset toward higher density at a given R a . 
This offset probably reflects the morphology-density relation: 
For a given Mg 2 , elliptical galaxies reside in denser environ- 
ments than lenticular galaxies. In addition, assuming that our 
sample of elliptical galaxies is morphologically homogeneous, 
ie. not contaminated by SO galaxies, the fact that they are also 
segregated indicates that the population-density relation may 
be more than a by-product of the morphology-density relation. 

Negative R a values may be due either to a low Mg 2 or 
symmetrically to a high <tq ■ While the first hypothesis is the 
most favored in the frame of the usual interpretation, the second 
may result from the dynamical evolution of galaxies harassed 
by gravitational encounters. In order to check this possibility 
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Fig. 1. Mean density of the environment log(p)as a function 
of the residuals R a to the Mg 2 - <to relation, R a . The filled 
symbols represent the whole sample, while the open symbols 
are the elliptical galaxies alone. Each point represents average 
log(p) for 80 (whole sample) or 40 (elliptical) galaxies. The 
error bars are the 1-er uncertainty on the position of each point. 



we will analyze the residuals to the FP: A low Mg 2 or a high 
do will result in opposite residuals. 

2.3. Residuals to the Fundamental Plane 

We performed a fit of the FP using the subsample of 291 bona- 
fide elliptical galaxies (17 galaxies were excluded from the 308 
galaxies sample because of the lack of reliable photometry): 



R f = 2 log (To + 0.2(1 + 2(3)M B +0.2^ e 



(2) 



Rf is the residual to the FP, o-q , the central velocity dispersion 
(in km sec -1 ) Mb and /i e respectively the asymptotic magni- 
tude and mean surface brightness within the effective aperture 
in the B band, (3 and r\ are the free parameters. Their best fit 
values are: (3 = 0.20 ± 0.02, r\ = -3.2 ± 0.2, similar to the 
values previously reported in PS. 

Note that we fitted the classical FP, ie., we did not include 
the additional terms accounting for the stellar population, the 
rotational support and the non-homology of the spatial struc- 
ture (see PS). Since we precisely want to use the residuals to 
the fundamental plane as a parameterisation of the stellar pop- 
ulation, we did not include the first scaling relation in the FP 
equation. Including the two latter would have dramatically re- 
duced the size of the sample and is not useful to the present 
goal. 

In Fig. H we present the average of log(/?) as a function of 
Rf. Each point on the figure represents 40 galaxies. There is a 
(somewhat marginal) tendency for the galaxies with the more 
negative residuals to be found in lower density environments. 

Splitting the sample in two parts, we find for Rf < 
-0.1, log(p) = -0.35 (154 galaxies), and for R f < -0.1, 
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Fig. 2. The mean relative density of the environment 
log(p) plotted as a function of the residuals Rf to the Funda- 
mental Plane. Each point corresponds to a bin of 40 galaxies, 
the error bars are the 1 — a uncertainty on the position of the 
points. 

log(p) = —0.19 (264 galaxies). The uncertainty on the mean 
log(p) being about 0.06, the difference is significant at the 3 a 
level. 

Repeating this on arbitrary extracted subsamples of half the 
total size produce the same result. Therefore, we conclude that 
the trend is significant. Since there is no reason why the mea- 
surement errors could be connected with the environment, we 
conclude that the relation is physically significant. 

Combining the results from the R a and Rf analyses, we 
confirm the existence of a population segregation, and extend 
it to the low-density environments. 

3. The origin of the population segregation 

Our result fits in the frame of the current paradigm suggest- 
ing a reduced rate of star formation in high density environ- 
ments. The morphology-density relation, the HI deficiency of 
cluster spirals or the relative isolation of shell galaxies are clas- 
sically interpreted as the result of past gravitational interactions 
which both stimulated the star formation earlier in the life of the 
galaxies and stripped the dense gas, prohibiting present-epoch 
star formation. 

We will now try to extract the fraction of the population 
segregation due to recent merging or accretion of a gas-rich 
companion, which cannot be accounted for by the morpholog- 
ical segregation. 

3.1. Evidence for recent merging events and star formation 
triggered by gravitational encounters 

Apart the violent cases of star forming galaxies often asso- 
ciated with strong interactions or mergers, the "weak" inter- 
actions have also been related with the presence of a young 



stellar component. In partic ular, Schweizer et al. (199C) and 
Schweizer & Seitzer ( 1992 > found clear correlations between 



related with the residuals to the FP (peculiar velocities in his 
terminology). To summarize, the presence of a "young" sub- 
population is clearly associated with merging events or grav- 
itational encounters. Hence, this population effect is expected 
to depend on the environment. 



Indeed, Guzman et al. (1992) found a difference between 



the galaxies located in the outer and inner regions of the 
Coma cluster. The former have negative residuals (i.e. they 
are younger or less metallic). Bower et al. (1992) found a dif- 



ference between the Coma and Virgo clusters, but this could 
also be an observational effect due to the difference in the pro- 
jected slit sizes. Recently J0rgensen & J0nch-S0rensen i 



showed that the colors and the absorption-line indices of E and 
SO galaxies belonging to the poor cluster S 639 indicate that 
the stellar populations in these galaxies are probably younger 
(or less metallic) than those in rich clusters. This is in agree- 
ment with the re sults by de Carvalh o and Djorgovski ( 1992 ), 
J0rgensen ( |l997| ) and Bernardi et al. ( |l998| ). 

It is important to note that these previous works searched 
for mean-age differences as a function of the environment. For 



instance, Bernardi et al. (1998) found a zero-point difference 
on the Mg 2 - <jq relation of <5(Mg 2 ) = 0.007 between cluster 
and field ellipticals, compatible with our comparison between 
the high- and low-density subsampl es. As suming a single burst 
stellar population model (Worthey 1994 ), this corresponds to 
an age difference of about 1 Gyr. 

In the present work, we found a much stronger difference 
when grouping the galaxies according the R a or Rf residuals. 
The galaxies containing a younger stellar population are clearly 
found in lowest-density environment, but, in any environment, 
the majority of the early-type galaxies have only an old pop- 
ulation. The rare galaxies containing yo ung st ars form the tail 
of negat ive re siduals noted in Burstein ( 1988 ) and Prugniel& 
Simien ( 1996| ). This same skewness is also apparent in the dis- 
tribution of the color-magnitude diagrams. It is also reminis- 
cent of the population of blue galaxies in cl usters (Butcher- 
Oemler effect; Aragon-Salamanca et al. 1991 ). If we compare 
the median of R a , less sensitive to the skewness, instead of the 
average, the difference between the high- and low-density sub- 
samples vanishs. 

We also tried to connect the residuals with the small-scale 
density of the environment by weighting the galaxies accord- 
ing to their separation on the sky when computing the density, 
with the scope of defining a "strength of the tidal field". We 
experimented different weightings, but failed to find evidence 
for an effect related with the "tidal" field. Only the large scale 
environment correlates with the stellar content. 

The population segregation may still be contaminated by 
the mo rpholo gical segregation. For e xample, Kuntschner & 
Davies ( 1998| ) and Mehlert et al. ( |1998| ) find the ellipticals to be 
coeval in Coma and Fornax, and detect the presence of residual 
star formation in lenticulars only. 



the anomalies on the colors, Mg 2 , H(3 and CN, and the fine 
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3.2. Subtracting the morphological segregation 

If the morphological segregation results from the morpholog- 
ical segregation, R a should also correlate with any other pa- 
rameter dependent on the morphological type. Such parameter 
could be the signature of a disk, as detected from image anal- 
ysis (Scorza and Bender 1995) or from the rotational support 
(Vmax/co , where V max is the observed maximum rotational 
velocity). Both of these parameters are affected by projection 
effects and are available only for a restricted subsample. We 
will use Kiax/c^o , taken from HYPERCAT. 

In our sample, Vmax/oois available for 272 galaxies: 157 
ellipticals and 115 lenticulars, the mean values are respectively: 
< V max /a >= 0.27 ±0.02 and 0.43 ±0.04 (The correspond- 
ing mean log(p) are: -0.46 ± 0.04 and -0.50 ± 0.04). The 
difference is clear and shows that V max /ao is fairly correlated 
with the morphology, but the morphology-density relation is 
here very marginal. 

Because of the smaller statistics the \og(p) - R a correlation 
is much noisier than in Fig. |lj but the effect is still marginally 
observed in the subsample: The mean log(p) for R a < is 
-0.50±0.03 (162 galaxies) and -0.44±0.04 for R a > 0(114 
galaxies). The corresponding mean Vmax/oo are: 0.328±0.023 
and 0.288 ± 0.025. The difference is marginally significant 
0.040 ±0.035 

The Vi nax / (To - log(p) correlation, represents the morphol- 
ogy - density relation, it is fitted as: log(p) — 0.25V max /a . 
It predicts a difference in log(p) between the two subsamples, 
R a < and R a > 0, of: 0.010 ± 0.015. This is about 1/4 
of the observed effect. However, taking into account the large 
uncertainties, it is not possible to reject the possibility that the 
totality of the stellar population segregation may be due to the 
morphology - density relation. 

Repeating the analysis on Rf r eprod uces the dynamical ef- 
fect found by Prugniel & Simien ( 1994 ) and does not help to 
understand the origin of the population segregation. 



4. Conclusion 

Analyzing the residuals to the Mg 2 - <7 and Fundamental 
Plane relations as a function of the large scale density of the 
environment, we find a segregation of the stellar population. 
The early-type galaxies with a line strength index Mg 2 weaker 
than expected from their value of gq are preferentially found 
in low-density regions. These galaxies are likely to contain an 
excess of young stars, presumably formed after a gravitational 
encounter or merging event that occurred in the last gigayears. 
This extends previous results to very sparse environments. 

However, we cannot rule out that this effect is a by-product 
of the morphological segregation. Indeed, even our subsam- 
ple of bona-fide ellipticals may still be contaminated by some 
lenticular galaxies seen almost face-on or with a weak disk, 
and the observed segregation could be due to residual star for- 
mation in these galaxies expected to be found mostly in low- 
density regions. We tried to use the observed rotational sup- 
port, Vmax/co , to parameterize the morphology, as it traces 



the presence of a disk. The morphology-density relation seems 
to account for part of the segregation observed on stellar pop- 
ulation, but the data are still too noisy to exclude that it is the 
totality. New data on a larger sample are needed. 
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